Neutrino Luminosity and Matter-Induced Modification of 
Collective Neutrino Flavor Oscillations in Supernovae 

John F. Cherry/'^ Meng-Ru Wu,^ J. Carlson,^' ^ Huaiyu 
Duan,^'^ George M. Fuller/'^ and Yong-Zhong Qian^ 

^Department of Physics, University of California, 
^ San Diego, La Jolla, California 92093, USA 

^ ^Neutrino Engineering Institute, New Mexico Consortium, 

Los Alamos, New Mexico 87545, USA 

D 

C/^ '^School of Physics and Astronomy, 

^ University of Minnesota, Minneapolis, MN 55455, USA 

I— I ^Theoretical Division, Los Alamos National Laboratory, 

JjH Los Alamos, New Mexico 87545, USA 

4^ ^Department of Physics and Astronomy, 

Q University of New Mexico, Albuquerque, New Mexico 87131, USA 

^ (Dated: September 27, 2011) 

' ' Abstract 

>■ We show that the bump in the electron number density profile at the base of the hydrogen 

in 

envelope in 0-Nc-Mg corc-coUapsc supernovae causes an interesting interplay between ncutrino- 
electron and neutrino-neutrino forward scattering effects in the flavor evolution of low-energy 

ON 

O in the neutronization burst. The bump allows a significant fraction of the low-energy to survive 

^ by rendering their flavor evolution nonadiabatic. Increasing the luminosity of the neutronization 

> 

burst shifts the bump-affected to lower energy with reduced survival probability. Similarly, 

X 

H lowering the luminosity shifts the bump- affected neutrinos to higher energies. While these low 

energy neutrinos lie near the edge of detectability, the population of bump-affected neutrinos has 
direct influence on the spectral swap formation in the neutrino signal at higher energies. 
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I. INTRODUCTION 



Stars of ~ 8-10 Mq (Mq being the mass of the sun) develop 0-Ne-Mg cores at the end 
of their evolution. Capture of electrons by the Ne and Mg isotopes initiates the gravi- 
tational collapse of the core, which eventually produces a supernova and leaves behind a 
neutron star [U |2]. These 0-Ne-Mg core-collapse supernovae are the only case for which 
the neutrino-driven mechanism has been demonstrated to work by different groups [3|, H] . 
The success of this mechanism is largely due to the steep fall-off of the matter density 
above the core. As shown by Refs. [3HH], this special density structure also provides a venue 
where neutrino flavor transformation occurs under the influence of both neutrino-electron 
and neutrino-neutrino forward scatterings [51422]. In particular, the neutronization burst, 
which consists predominantly of Ve and signifies the breakthrough of the neutrino sphere by 
the supernova shock, experiences interesting flavor evolution including collective oscillations 
for the neutrino flavor mixing parameters found by experiments [23] . 




In this paper we explore another special feature of the matter structure in 0-Ne-Mg 
core-collapse supernovae in connection with flavor evolution of the neutronization neutrino 
burst. The hydrogen envelope has an electron fraction of Ye ~ 0.85. In contrast, the material 
below the envelope has ~ 0.5, reflecting weak interaction-induced neutronization during 
pre-supernova evolution. As the matter density p is a continuous function of radius, this 
produces a bump in the electron number density rig = pY^N^ {N^ being Avogadro's number) 
at the base of the hydrogen envelope, shown explicitly in Ref. 0. As noted in Refs. [5HH], for 
the normal neutrino mass hierarchy, this bump renders flavor evolution of the low-energy Ue 
in the neutronization burst nonadiabatic, giving rise to substantial survival probabilities for 
these Ug. Here we show that this bump facilitates an interesting interplay between neutrino- 
electron and neutrino-neutrino forward scattering in the flavor evolution of the low-energy 
Ug, and we show how this influences the collective oscillations of neutrinos at higher energies. 

II. NEUTRONIZATION BURST NEUTRINOS 

We assume a pure z/g burst emitted from the neutrino sphere at = 60 km with a total 
luminosity L^, = 10^^ — lO^^ergs"^ and a normalized spectrum 
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where we take r],y = 3 and T^, = 2.75 MeV [corresponding to an average z/g energy {E,y) = 
Fz{rjy)Ty / F2{rjy) = 11 MeV at emission]. Here 



In the single-angle approximation, the neutrino-neutrino forward scattering potential can 
be written in terms of an effective total neutrino number density aX r > R^, 



where the approximate equality holds for r ^ R^. 

For the purposes of this study we have chosen the following neutrino mixing parameters: 
neutrino mass squared differences Arrig = 7.6 x 10~^eV^ and Aml^j^ = 2.4 x 10~^eV^; 
vacuum mixing angles 6*12 = 0.59, 6*23 = vr/4, ^^13 = 0.1; and CP-violating phase 6 = 0. We 
have also chosen to use the single-angle approximation for these calculations, where neutrinos 
emitted at different angles relative to the surface of the neutrino sphere are assumed to 
have the same flavor evolution history as neutrinos from a single, representative emission 
trajectory. 

Here we will concentrate on the normal neutrino mass hierarchy, because previous work 
has shown that mixing at the Aml^^ scale with this heirarchy produces interesting collective 
neutrino flavor oscillations in 0-Ne-Mg core-collapse supernovae [5H8] . For the atmospheric 
neutrino mass doublet with an inverted hierarchy, neutrinos in the neutronization burst do 
not experience any flavor transformation. 

The results from single-angle simulations shown in Figures [T] and [2] demonstrate that 
the He proflle with the bump gives rise to a very clear flavor transformation signature. 
Depending on the luminosity, a population of Ue below 8 MeV have large probabilities to 
transform between neutrino mass states (have large hopping probabilities). For example, the 
medium range luminosity case, = 8.0 x lO^^ergs"^, exhibits a peak hopping probability 
for the bump-affected neutrinos of ~ 80% at = 4.5 MeV, shown in Figure [2j In the 
extreme case of L„ = lO^^ergs"^, the hopping probability is ~ 15% for E,, = 0.5 MeV, 
although this is hard to see in Figure [l|(but see Figure 2h. 
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FIG. 1: The final neutrino mass state emission energy spectra for calculations of the flavor trans- 
formation in the neutronization neutrino burst of an 0-Ne-Mg core-collapse supernova. Each panel 
shows the results for a different possible burst luminosity, ranging from = 10^^ - 10^^ 

erg s , 

with identical Fermi-Dirac energy distributions. 

III. ANALYSIS OF FLAVOR EVOLUTION OF THE s 

As discussed in Refs. [5H8] the flavor evolution of the Vf. flux in this case is governed by 
'^"^atm ^13- Although the numerical calculations we have conducted employ full 3x3 
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FIG. 2: The probability for electron neutrinos in the neutronization neutrino burst of an 0-Ne-Mg 
core-collapse supernova to hop out of the (initial) heavy mass eigenstate, 1 — Ph, plotted as a 
function of inverse neutrino energy. Each panel shows the results for a different possible burst 
luminosity, ranging from L^, = 10^^ — 10^^ ergs~^, with identical Fermi-Dirac energy distributions. 

flavor mixing, in the following analysis, we focus on the 2-flavor mixing for these z/g with 
= 5m^^jjj and 9^ = 613 ■ Further, we adopt the single-angle approximation, as this has 
been shown to be surprisingly effective in providing qualitative understanding of the results 
from multiangle simulations. 
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Using the notation introduced in Ref. [23], we can represent a Ve of energy Ehy a. neutrino 
flavor isospin (NFIS) with uj = 5m? /2E. The evolution of is governed by 

d 



dr 



poo 

wHv + He - At(r) / s^>fy{E^,)dE, 
Jo 



(4) 



where Hv = cos2^ve2 — sin2^ve^, He = ~V^GFneir)el, yu(r) = 2A/2G'i?ny(r), and E^' = 
5m? /2uj' . Here e^. and e\ are the unit vectors in the x and z directions, respectively, of the 
neutrino flavor space. For convenience, we deflne 

9{^) = -^UE^) (5) 



and 



It follows that 



/»0O /'OO 

/ s^fy{E^)dE^ = / s^g{uj)du. (6) 
Jo Jo 

POD 

-S = / ug{u)s^du X Hv + S X He. (7) 
Jo 



d 
dr 

As g{u) is concentrated in a flnite range of cu, to zeroth order we approximate g{uj) ~ 
5{co — (w)), where (u) = ujg{uj)du is calculated from the actual g{uj) in Eq. (5). Then 
the zeroth-order mean fleld S'-''^ can be obtained from 

|;S(°) = S(°) X [(u;)H. + He] = S(°) x Hmsw- (8) 

The evolution of S^°^ is the same as that of a Ve with -Emsw = 5m'^ /2{uj) = 8.53 MeV 

undergoing the usual MSW effect. With this, we can approximately solve the evolution of 
by employing 

^s, ^ X [ioU, + He - /x(r)S(°)] . (9) 

As the heavy mass eigenstate essentially coincides with at high densities, but the light 
mass eigenstate is predominantly z/e at low densities, the survival probability of an initial 
Ue is approximately 1 — Ph, where Ph is the probability for remaining in the heavy mass 
eigenstate. 

A. Dependence on L,y 

We can go further by using the zeroth-order mean fleld S^"-* to understand how the flavor 
evolution of the low-energy z/g depends on L^. As discussed below. Equations [8] and |9] imply 



that neutrinos with u ^ (w) will experience an MSW resonance before the resonance of 
S*^*'-*. These higher frequency neutrinos may pass through multiple resonances created by 
the matter potential bump. 

Based on the MSW effect, S^°^ corresponding to -Emsw = 8-53 MeV goes through the 
resonance after the low-energy u^- Assuming adiabatic evolution of S'-'^-' before the resonance, 
we can take 

Hmsw ^ 



where 



S(o) 



cos 29r^ 



sin 29r. 



2|H 



MSW I 



--(cos26l„e{ - sin26'„e^). 



(cj) cos2^^v — v^G^rie 

((w) cos 26^ - ^GpUef + ((w) sin2^v)^ 
{(jj) sin 26'v 



\{uj) cos 29^ - ^GpUeY + ((w) sin 20V )^ 
The evolution of before the resonance of S*^°^ is then governed by 
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dr 



s,,. X 



[ijj cos 29^ — v^G^ne + ^ cos 29m)^{. — {u sin 29^ + ^ sin 29m)e 



fx 



The above equation shows that s^^ goes through the resonance when 



u cos 29, = V2G 



COS 29r, 



IH, 



B. 



(13) 
(14) 

(15) 



Note that cos20m < before the resonance of S'^'^^ [see Eq. (11)] and therefore B > 0. 
Consequently, the energy of those {E^ = 5rn? /2uj) that go through the resonance at the 
bump in the rif. profile decreases as Ly, and hence B, increases. This trend can be seen in 
every frame in Figures [T] and |2j As the neutrino luminosity is increased, the peak energy of 
the population of low energy neutrinos that hop out of the heavy neutrino mass eigenstate 
as a result of the bump decreases. 

We can also qualitatively understand why increasing produces a decreasing survival 
probability of the bump-affected z/g. The Landau-Zener probability for hopping from the 
heavy to the light mass eigenstate after the resonance is 

vr 5im? sin^ 29^ 



where 



Phop = exp 



4 E cos 29. 



rfln(|He| +5) 
dr 



(16) 



(17) 



is the scale height of the total flavor-evolution potential at the resonance position. Crudely 
we have 1 — Ph ~ -Phop- As increasing Lj^ shifts the resonance energy window to lower Ey at 
the bump in the rie proflle, Phop decreases because flavor evolution through the resonance 
tends to be more adiabatic for lower-energy neutrinos [see Eq. (16)]. In addition, as B 
decreases much more slowly than |He| with radius, l-Lres becomes larger when the contribution 
from B increases with Ly. This also reduces Phop [see Eq. (16)]. 

Furthermore, equations [8] and [9] also imply that a sufficiently large neutrino-neutrino 
scattering potential will cause neutrinos with oscillation frequencies roughly equal to or less 
than (u) to follow the evolution of S^"-* as this vector moves through resonance. To illustrate 
this, we choose to define the angle a as the angle between Hmsw and either S*^°^ or S. The 
heavy mass eigenstate survival probability Ph of the collective ensemble of neutrinos that 
follow the evolution of S^"-* is related to a by 

Ph = 1-^(1 + cosa). (18) 

Figure |3] shows the evolution of a for the simulations with relatively high luminosities. From 
the figure it can be seen that the evolution of S for these luminosities is qualitatively similar 
to that of S(°). 

Interestingly, the final alignment angle, a, for the collective neutrino isospin vectors is 
slightly larger than it is for S^°^ This means that the collective NFIS's are more closely 
aligned with — Hv than 8*^°-* is. The highest luminosity simulation, with Ly = lO^^ergs"^, 
has the collective NFIS that is most closely aligned with S'-^-', and the reason for this can 
be found in Eq. M In the limit of (r) ^ |He| ^ a; the individual will orbit exclusively 



around S*^''^ and follow it through resonance. However, it can be seen from Eqs. 11 and 13 that 
only neutrinos with u = (u) go through resonance at the exact position where cos26'jn = 0. 
Neutrinos following the evolution of S^"-* will still experience some fraction of the neutrino 
self-coupling potential, although at resonance |He| ^ B for neutrinos that track 8*^*^^ This 
results in a small increase in T^resj which slightly lowers the overall hopping probability and 
slightly increases a. 

When the neutrino luminosity is moderately lower, the same basic phenomenology is 
observed. Figure |4] shows the evolution of a for the simulations with moderate luminosities, 
0.6 — 1.0 X lO^^ergs"^. The collective NFIS 8 for these simulations still tracks roughly the 
evolution of 8^°\ although deviations become more pronounced as the neutrino luminosity 
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FIG. 3: High luminosity evolution (Lq = 10^'^ ergs ^): The opening angle a between the collec- 
tive NFIS S(o) and Hmsw, plotted function of | Hmsw I / |Hv| as the system moves through 
resonance. The idealized NFIS (solid line) shows the evolution of S^'^^ in the ideal, strong neutrino 
self-coupling case. The dashed line, dot-dashed line, and dotted line show the evolution of S as 
calculated for neutrino luminosities IOLq, \/10^05 Lo respectively. 
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FIG. 4: Moderate luminosity evolution (Lq = 10^^ erg s^^): The opening angle a between the 
collective NFIS S*^*'^ and Hmsw, plotted as a function of |Hmsw|/|Hv| as the system moves 
through resonance. The idealized NFIS (solid line) shows the evolution of S^*^-* in the ideal, strong 
neutrino self-coupling case. The dashed line, dot-dashed line, and dotted line show the evolution 
of S as calculated for neutrino luminosities Lqi 0.8 Lq, and 0.6 Lq respectively. 
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decreases. Counter- intuitively, the final alignment of the lower luminosity S's is closer to 
that of S^"-* than in the calculations with = 10^^ erg s^^. This effect originates in the 
contribution of the bump-affected neutrinos to the integrals in Eqs. |6] and [7j From Figure [T] 
one can see that the population of bump-affected neutrinos has grown appreciably in this 
luminosity range, comprising 7 — 10 % of all neutrinos. These bump affected neutrinos are 
not connected in a coherent fashion to the flavor evolution of S, but they are predominantly 
aligned with the -t-H^ axis (they are predominantly z/2). This means that they will tend to 
drag the alignment of S closer to the -|-Hv axis, which systematically moves the final value 
of a lower. 

Opening Angle a Relative to Hmsw 

3.5 
3.0 
? 2.5 

2.0 
1.5 
1.0 

10"' 10" 10' 10^ 10^ 10" 

|Hmsw|/|Hv| 

FIG. 5: Low luminosity evolution (Lq = 10^'^ ergs^^): The opening angle a between tlie collective 
NFIS S(°) and Hmsw, plotted as a function of | Hmsw I / |Hv| as the system moves through res- 
onance. The idealized NFIS (solid line) shows the evolution of S*^^^ in the ideal, strong neutrino 
self-coupling case. The dashed line, dot-dashed line, and dotted line show the evolution of S as 
calculated for neutrino luminosities 0.4 Lq, 10~^/^Lo, and 10^^ Lq respectively. 

For luminosities below L,y = 6.0 x 10^^ erg s~^, the magnitude of the neutrino self-coupling 
potential drops below (w) = 1.4 x 10~^^MeV prior to reaching the resonance region of 8^*^^ 
This means that many neutrino states will undergo an MSW-like flavor transformation prior 
to reaching this region of the supernova envelope. As a result, these neutrinos, including 
those with u = (u) , will not participate in the collective flavor oscillations we have described. 
In turn, this means that the approximation g (u) ~ 5 (w — (u)) is not justified in this case. 
Ultimately this approximation breaks down because the evolution of the the neutrino state 

10 




with UJ = (u) is not coherently related to the evolution of other neutrino flavor states. 

The progressive breakdown of this approximation with decreasing neutrino luminosity 
can be seen in Figure |5} The motion of the vector S for each calculation deviates widely 
from the motion of S^^\ Furthermore, the flnal alignment angle a for each S is no longer 
related to the actual hopping probability for neutrinos in those calculations. The hopping 



probability inferred from Figure |5] and Equation 18 differs dramatically from the actual 
hopping probability observed in the calculations shown in Figures [T] and [2] for low L,^. The 
discrepancies are APhop = 0.16, 0.35, 0.48 for the calculations with = 4.0, a/TO, 1.0 x 
]^Q52gj,gg-i respectively. 



B. Spectral Swap Formation 

If the luminosity is large enough, it can be seen that after the neutrinos in these calcu- 
lations have passed the resonance region, the neutrino self-coupling potential becomes the 
dominant term in the neutrino forward scattering potential. Neutrinos which have u < ^{r) 
fall into a form of collective flavor oscillations known as the Regular Precession mode. The 
Regular Precession mode is typifled by the collective ensemble of neutrinos rotating with a 
common frequency, about the axis of the vacuum mass basis, 

= Up, (s^ X Hv) . (19) 

This collective oscillation has the feature that it conserves an effective lepton number (or 
"energy") of the ensemble of neutrinos. While this lepton number has a more complicated 
general expression, in the particular case of the neutronization neutrino burst where the 
initial flux of neutrinos is nearly pure z/e, the conserved lepton number is simply oc n^Pn for 
neutrino mixing at the atmospheric mass scale. 

In the initial stages of neutrino flavor transformation, this lepton number is not conserved. 
However, thereafter the Regular Precession mode flxes the total number of neutrinos in mass 
state 3. This gives the criterion for the precession frequency, Up,, for the system, 

/•oo rco 

/ Ph (to) g{uj)duj= g{uj){PH (a;) - [P2 {uj) + Pi (c^)] }duj, (20) 

J Wpr J 

where Pi {u) , P2 (w) are the probabilities of a neutrino with oscillation frequency u to be 
in the instantaneous mass eigenstate 1, or 2 respectively. 
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It is this precession frequency that sets the energy of the spectral swap (in this case 
between mass state 3 and mass state 2). As the magnitude of the self coupling drops, 
neutrinos with oscillation frequencies in the range fi{r) > u > Upj- will participate in the the 
Regular Precession mode and will align with mass state 3, while neutrinos with u < will 
be aligned with mass state 2. The final results of this process can be seen in Figure [Tj For all 
but the least luminous calculation a spectral swap forms, with -Eswap = Sm'^/2upr, close to 

~ 15 MeV. The precise location of -Eswap depends on the details of flavor transformation 
due to the motion of S'-'^-' and the bump affected neutrinos. Broadly speaking, a smaller 



Ph found from Eq. 18 will lower the swap energy of the final neutrino energy spectra by 



reducing the value of the integral on the right side of Eq. [20| However, a larger population 
of bump affected neutrinos will move the swap energy to higher values (smaller Up^) by 
reducing Ph (i^) 9 for large u. 

It is important to note that the spectral swap between mass states 3 and 2 can sometimes 
form even when the coherent flavor evolution of neutrinos has broken down deeper in the en- 
velope. As discussed in the previous section, for the calculations with < 6.0 x 10^^ erg s~^ 
collective neutrino flavor transformation breaks down in the resonance region because 
/i (r) < {cu). However, from Figure [T] it can be seen that a mass state 3/2 swap still forms 
successfully for = 4.0 x lO^^ergs"^ and Li, = vTO x lO^^ergs"^. Swaps between mass 
state 3 and 2 form for these two models where the luminosity is low because the neutrino 
self-coupling is still large compared to u for high energy neutrinos, specifically fi (r) > 
after the resonance region. This allows the high energy neutrinos to briefly form a Regular 
Precession mode before fi (r) decreases further with radius and flavor transformation in the 
Sml^^ mixing sector stops. 

The swaps are incomplete for the calculations with L^, = 4.0 x lO^^ergs"^ and Li, = 
\/lO X lO^^ergs"^. The calculations for these cases show small populations of 1^2 neutrinos 
below -Eswap which are not bump affected. They also show small populations of neutrinos 
in z/3 with energies above -Eswap- This phenomenon arises because of the extremely short 
lifetime of the collective precession at these luminosities. We define rstop to be the distance 
between the end of the resonance region and the point at which fi (r) = Wpj.. Comparing 
that to the oscillation length of neutrinos in the regular precession mode, /qsc = Syr/wpr, 
we find that for these two cases rstop ~ ^osc- Clearly, a spectral swap cannot fully form if 
the Regular Precession mode ceases before it can complete several full oscillations. It is 
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interesting, however, that the swaps in these two calculations are as pronounced as they are 
given the rapid truncation of the collective neutrino oscillations. 

For the lowest neutrino luminosity, = 1.0 x lO^^ergs"^, there is no spectral swap 
observed in Figure [l] between mass states 3 and 2. In this case, /i (r) < Upr even before the 
system finishes MSW-like flavor transformation. No collective oscillation can proceed in the 
^"^atm mixing sector for this case. 

While we have focused entirely on the (5m^^j^ mixing sector in the sections above, it 
should be pointed out that the StUq mixing sector is completely indifferent to the range of 
luminosities that we have explored. The neutrino flavor mixing energy scale for the SitLq mass 
state splitting is ~ 30 times smaller than that of the atmospheric mass state splitting. The 
model of neutrino flavor transformation outlined above is quite robust for the solar mixing 
sector, with fi (r) > {uj)q and fi (r) > (wpr)Q for all of the neutrino luminosities that we 
consider. The spectral swap between mass states 2 and 1 is created by the Regular Precession 
mode in this mixing sector. While -EswapQ varies greatly for the different calculations in 
Figure [1} this swap energy is only changed by variations in flavor transformation in the 
^"^atm sector. The ratio of 1^2 /i^i neutrinos is identical for all of the calculations shown in 
Figure [l} A curious consequence of this is that the spectral swap energies move closer and 
closer together as the luminosity of the neutronization burst decreases, until ultimately the 
swap between mass state 3 and 2 disappears altogether. This behavior is evident in Figure [1} 



IV. DISCUSSION AND CONCLUSIONS 



We have shown that the bump in the electron number density proflle at the base of 
the hydrogen envelope in 0-Ne-Mg core-collapse supernovae causes an interesting interplay 
between neutrino-electron and neutrino-neutrino forward scattering effects in the flavor evo- 
lution of low-energy z/g during the neutronization burst epoch. The bump allows a signiflcant 
fraction of the low-energy z/g to remain in the electron flavor state, i.e., enhancing their sur- 
vival probability. It does this by rendering their flavor evolution nonadiabatic. Additionally, 
we have found that increasing the luminosity L^, of the neutronization burst shifts the bump- 
affected z/e to lower energy with consequently reduced survival probability. Finally, we have 
found that the flavor states of the bump affected low energy neutrinos impact the spectral 
swap forming behavior of the collective oscillations later on. This opens up the possiblity 
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that the presence of bump affected low-energy v^. of ~ 3-5 MeVcan be detected through the 
spacing of the spectral swaps. However, this may not be necessary as some neutronization 
burst signals may produce bump affected neutrinos at energies that are accessible by Earth 
based detectors. This may also prove to be an interesting secondary probe of the burst 
luminosity L^, of an 0-Ne-Mg core-collapse supernova. 

While the zeroth-order mean field proves to be rather useful in understanding the flavor 
evolution of the low-energy v^.^ it is clearly inadequate in providing a quantitative description 
of the flavor evolution if the v^. burst luminosities are on the low side of the expected emission 
(see Figure |5]). In particular, it cannot provide a good estimate for the energy at which the 
spectral swap occurs in the low luminosity limit. We note that the spectral swap occurs 
rather robustly at a fixed energy once Ly exceeds ~ vTO x 10^^ erg/s (see Figure 
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